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Abstract. We are concerned with viscous profiles (travelling waves and steady 
solutions) for mixed hyperbolic-parabolic systems in one space variable. For 
a class of systems including the compressible Navier Stokes equation, these 
profiles satisfy a singular ordinary differential equation in the form 

<»•" %=m Fm ' 

Here U takes values in R d and F : — > M d is a regular function. The real 
valued function C(£-0 is as we H regular, but the equation is singular because 
C(U) can attain the value 0. We focus on a small enough neighbourhood 
of a point U satisfying F(U) = 0, C(^0 = 0- From the point of view of the 
applications to the study of hyperbolic-parabolic systems this means restricting 
to systems with small total variation. 

We discuss how to extend the notions of center manifold and of uniformly 
stable manifold. Also, we give conditions ensuring that if £(U) ^ at t = 
then C(t/) ^ at every t. We provide an example showing that if f (U) becomes 
zero in finite time then in general the solution U of equation I IO.U is not 
continuously diffcrentiable. 



1. Viscous profiles for mixed hyperbolic-parabolic systems 

We are concerned with systems in the form 
(1.1) E(u)u t + A(u, u x )u x = B(u)u xx . 

Here the function u takes values in M. N and depends on two scalar variables, t and 
x. The matrices E, A and B have all dimension N x N. The conservative case 



u t + f(u) x = (B(u)u, 



is, in particular, included in the formulation (II. ip . 

In the following, we will focus on mixed hyperbolic-parabolic systems, i.e. we 
will assume that matrix B in (jl.ip is singular. This case case is interesting from 
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the point of view of applications. Indeed, it is the case of the compressible Navier 
Stokes equation in one space variable: 



Pt + (pv) x = 
(1.2) J W + ^+p)^^^ 



1 ., 
-pv +pe + p 



{^k6 x + vvv x ^j 



Here the unknowns are p(t, x), v(t, x) and e(t, x): the function p is the density of 
the fluid, v represents the velocity of the particles in the fluid and e is the internal 
energy. The function p = p(p, e) > is the pressure and satisfies p p > 0, while 9 is 
the absolute temperature and in the case of a polytropic gas satisfies 

e _ e(7~l) 
R ' 

where R is the universal gas constant and 7 is a constant specific of the gas. Fi- 
nally, v{p) > and k(p) > are the viscosity and the heat conduction coefficients 
respectively. 

In 12 , Kawashima and Shizuta introduced a set of hypotheses that are satisfied 
by the equations of the hydrodynamics and of the magnetohydrodynamics and that 
are frequently exploited to study the hyperbolic-parabolic system (jl.ip . In the 
following we will assume that the Kawashima Shizuta conditions are satisfied but, 
since we do not exploit them explicitly, we do not recall them. 

We are concerned with special classes of solutions of (jl.ll) . namely travelling 
waves and steady solutions. Travelling waves satisfy the ordinary differential equa- 
tion 

(1.3) [A(U, U') - aE{U)] U' = B(U)U", 
while steady solutions are solutions of the ODE 

(1.4) A(U, U')U' = B{U)U". 

In (| 1 . 3[) , a is a real constant which is usually called the speed of the wave. Note that 
from a solution U(y) of (|1.3p we can obtain a solution of the original hyperbolic- 
parabolic system (II. ip by setting u(t, x) = U(x — at). Moreover, any solution U(x) 
of (|1.4p is a steady solution of (jl.ip . i.e. a solution that does not depend on time. 
Also, in the following we will focus on the case <r, the speed of the travelling wave 
(|1.3p . is close to an eigenvalue of the matrix A(U, 0). Since in general is not an 
eigenvalue of A(U, 0), we keep the cases (jl.3p and (|1.4|) separated. 

It is known that the study of travelling waves and steady solutions can provide 
useful information to study the limit e — > + of the family of functions u e satisfying 

(1.5) E(u e )u e t + A(u e , eu%)u £ x = sB(u £ )u xx . 

The literature concerning this issue is very wide, so we just refer to Benzoni-Gavage, 
Rousset, Serre and Zumbrun 2 , to Rousset |14j . to Zumbrun [17] . and to the rich 
bibliography contained therein. For a more general introduction to the parabolic 
approximation of hyperbolic problems we refer instead to the books by Dafermos 
[3 and by Serre [15] and to the references therein. 

Note that, when system (jl.ip is the Navier Stokes equation, the system we 
obtain formally setting e — in (|1.5p is the Euler equation. 
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In studying the viscous profiles of the Navier Stokes equation (11. 2| one encoun- 
ters a singular ordinary differential equation in the form 

M § - 

Here the unknown [7 is vector valued and has the same dimension as the function 
F. The function £ is real valued and the singularity of the equation comes from the 
fact that £ can attain the value 0. The link between (|1.2[) and (| 1 . 6[) is the following. 
Assume that we want to study the steady solutions of the Navier Stokes equation: 
we thus focus on (11.41). Set 



(1.7) 



P.r 



\Vx, 



After some computations (see Bianchini and Spinolo [4] for the details), we get that 
the steady solutions of the Navier Stokes equation satisfy 





(1.8) 



A 



21 



A22 



W 

z 



b(u) 



Here A22 and b are 2x2 matrices, an is a real valued function, the function A21 
takes values into K 2 and A t 21 denotes its transpose. The exact expression of these 



terms is not important here: 
Equation (|1.8[) gives 



we just point out that the matrix b is invertible. 



anvw 
A 2 iw - 



YA\ x z = 
A22Z = bz x 



Assume anv ^ 0, then we get 



(1.9) 



At 



anv 



= b 



-1 



.4 



A21A 



22 



21^21 



Note that the previous expression is well defined since the matrix b is invertible. 
Summing up, we get that the steady solutions of the Navier Stokes equation can 

z)\ ({U) = v and 



be written in the form (|1.6p provided that U = (p, v, e 

A\ x zja\\ 

F(U) 

A 22 v - A^iAlJaix 



\ 



J 



Some remarks are here in order: first, the exact expression of the function an is 
Pp/p 2 , where p p > is the partial derivative of the pressure with respect to the 
variable p. We can restrict to the case p is strictly positive and bounded away from 
zero: this implies that vacuum states are not assumed. We then have that the 
function an is well defined and does not attain the value zero. On the other side, 
v represents the velocity of the fluid and in general it can attain the value zero, 
which is the singular value for the equation satisfied by the steady solutions. 

Moreover, we underline that so far we have considered only steady solutions. 
However, also the equation of travelling wave profiles (| 1 - 3[) may become singular, 
in the sense of (II. 6|) . 

Also, the considerations carried on so far can be extended to a larger class of 
mixed hyperbolic-parabolic systems that do not satisfy a condition of block linear 
degeneracy defined in Bianchini and Spinolo [5]. 
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Remark 1.1. The reason why we have the term an in (|1.8|) is the following: 
instead of working directly on (|1.2j) , we consider an equivalent symmetric system 
in the form (see again Bianchini and Spinolo [4] for the details). System (|1.1|) 
is symmetric if 

A(u, 0)* = A(n,0). 

2. Invariant manifolds for a class of singular ODEs 

In the following we will thus focus on the singular ordinary differential equation 

<"> f " ck Fiir >- 

We will be concerned with the solutions U belonging to a small enough neighbour- 
hood of a value U satisfying C(t7) = 0, F(U) = 0. From the point of view of the 
applications to the analysis of system this implies that we are restricting to 

steady solutions and travelling waves having small enough total variation. Without 
loss of generality, in the following we will assume that U = 0. 

What we are interested in is the existence of locally invariant manifolds for 
(|2.1[) . If U takes values in M. d , a locally invariant manifold M. for (|2.1[) is contained 
in R d and satisfies the following: if uo 6 M, then the solution of the Cauchy 
problem 

1 :F(U) 




-- u 

belongs to M. is \t\ is small. In particular, we are interested in extending the notions 
of uniformly stable and of center manifold to the case of the singular ODE (|2.1[) . 
We recall here that a center manifold for the non singular ODE 

(2.2) ^ = G(U) U e K d 

at 

is defined in a neighbourhood of an equilibrium point U. Loosely speaking, a center 
manifold contains the orbits of (|2.2p that are globally bounded for t — > ±oo , more 
precisely \U(t)\ < S for every t. The constant S is the size of the neighbourhood 
and depends on the function G and on the point U. We refer to the book by Katok 
and Hasselblatt [11] for a complete discussion. Also, a presentation of the the most 
important properties of the notion of center manifold is given in Bressan's notes [6]. 
For the applications of the notion of center manifold to the study of the parabolic 
approximation of hyperbolic problems see for example Bianchini and Bressan [3] . 

The stable manifold of (|2.2|) contains the orbits that when t — > +oo converge 
exponentially fast to the equilibrium point U. We refer to the book by Perko [13] 
for a complete discussion. The notion of uniformly stable manifold is an extension 
of the notion of stable manifold. Assume that E is a manifold of equilibria for 
system (|2.2I) and assume that U <E E. The uniformly stable manifold relative to E 
contains the orbits of system (|2.2[) that for t — ► +oo decay eponentially fast to an 
equilibrium point in E. It is defined in a small enough neighbourhood of U and 
contains the stable manifold. However, in general the inclusion is strict since the 
orbits on the stable manifold can converge only to U, while when we are on the 
uniformly stable manifold the limit can vary on E. The uniformly stable manifold is 
sometimes called the slaving manifold relative to E. Its existence can be viewed as 
a consequence of Hadamard-Perron Theorem: we refer again to the book by Katok 
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and Hasselblatt [11] for the Hadamard-Perron Theorem, while a specific discussion 
on the uniformly stable manifold and its applications to the study of the parabolic 
approximation of hyperbolic systems can be found in Ancona and Bianchini [1] . 

It is known that, from the point of view of the applications to the study of the 
parabolic approximation (|1.5p . if we restrict to systems having small enough total 
variation it is interesting to focus on travelling waves lying on a center manifold 
and steady solutions lying on a stable or on a uniformly stable manifold. 

There is a rich literature concerning the family of systems 

(2.3) ^-=-F(U,e) (7el d . 

at e 

Here we just refer to the notes by Jones [10] and to the rich bibliography contained 
therein. In particular, |10j provides a nice overview of some papers by Fenichel ([9, 
[8] for example), whose ideas are exploited in the following. In (]2.3[) the singularity 
e is a parameter, e — > + . The main novelty here is that in (|2.1| we consider the 
case C(U) is a nontrivial function of the solution U itself. As a consequence, we 
have to take into account the possibility that ((U) ^ when t = 0, but £([/) attains 
the value in finite time. This may lead to a loss of regularity in the solution U. 
Consider, for example, the following system: 

dui/dt = —uiju\ 

d,U2/dt = —U2- 

It can be written in form (|2 . 1 1) provided that U = (tti, u 2 )*, ((U) = u\ and 

F{U) = ( ~ U2 



(2.4) 



The solution of |2l} is 
(2.5) 



t*i(*) = V /ui(0)+u 2 (0)(e 



«a(*) =«2(0)e-* 

Choosing u 2 (0) > Ui(0) > 0, one has that ((U) = ui(t) can attain the singular 
value for a finite t. Note that at that point t the first derivative du\j dt blows up: 
thus, in particular, the solution (|2.5[) of (|2.4[) is not C 1 . 

In the following, we will look for conditions that rule such a loss of regularity: 
this sounds reasonable in view of the applications to the analysis of the viscous 
profiles. Indeed, when considering the parabolic approximation (|1.5|) it seems rea- 
sonable to look for regular solutions. To prevent losses of regularity like the one in 
(|2.5p we need to be sure that if a solution U of (|2.1[) satisfies ((U) at t = 0, 
then C(U) 7^ for every t. In the following, we will state conditions ensuring that 
this property holds. 

Remark 2.1. To simplify the notations, in the following we will always assume 
that ((U) > at t = 0. The case ((U) < does not involve additional difficulties 
and can be tackled with techniques similar to those discussed here. 

3. Main results 

Summing up, our goals are as follows: we want to define locally invariant 
manifolds that extend the definition of center and uniformly stable manifold and 
we want to make sure that losses of regularity like the one in (|2.5p are ruled out. 
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We state our main results here, postponing to Section |4] the statement and the 
discussion of the hypotheses. We refer to [4] for a proof of Theorems 13.11 and 13.21 
Concerning the definition of a center manifold, we have the following. 

Theorem 3.1. Let all the hypotheses introduced in Section^hold. Also, assume 

that 

F(0) = C(0) = 0. 

Then there exists a manifold A4 C , defined in a small enough neighbourhood ofO and 
satisfying the following properties. The manifold M c is locally invariant for 

It contains all the solutions U satisfying \U{t)\ < 6 for every t. Here 6 > is 
a constant depending on the system. Also, if U(t) is an orbit lying on Ai c and 
C(U) ^ at t = 0, then ((U) for every t. 

Before discussing the extension of the notion of uniformly stable manifold, some 
remarks are in order. One can show (see Bianchini and Spinolo [4] for the details) 
that the solutions of (|3.1[) that lie on the manifold M c satisfy an equation which 
actually has no singularity in it. The situation is completely different if we consider 
the uniformly stable manifold. To see this, let us consider a trivial example: 

!du\jdt — —hu\ 
du 2 /dt = -u 2 /e 
de/dt = 0. 

Note that the subspace E = {(0, 0, e) : e £ K} is entirely made by equilibria of 
(I3.2p . If e > the uniformly stable manifold relative to E is the whole space R 3 , 
since any solution of (|3.2j) decays exponentially fast to a point in E. However, the 
first component goes like e -5 *: thus, the speed of exponential decay is bounded in 
e and u\ is not affected by the presence of the singularity. Conversely, the second 
component decays like e _t / e and hence the speed of exponential decay gets faster 
and faster as e — > + . The second component of the solution is thus strongly 
affected by the presence of the singularity and can be regarded as a fast dynamic, 
while the first component is a slow dynamic. Summing up, in (|3.2| any orbit lying 
on the uniformly stable manifold relative to E decomposes as the sum of a slow 
and a fast dynamic. 

This behaviour is somehow inherited by general non linear systems, in the 
following sense. If we perform the change of variable r = t/e, system (|3.2[) becomes 

dui/dr = — 5ux£ 

du 2 jdT = —U2 

de/dr = 0. 

We can thus single out the fast dynamics by saying that they are solutions of (|3.2p 
that are exponentially decaying with respect to the variable r. In the general non 
linear case, we consider system 



(3.3) 
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which is formally obtained from system (|3.1|) through the change of variables 
t = r(t) defined by 

' cZr 1 
(3 .4) J dt ~ Wit)} 

k r(0)=0. 

The /asi dynamics of (I3.1|l arc then the solutions of (13. 3j) that are exponentially 
decaying to zero with respect to the variable r. Note that, a priori, the change 
of variable Ij3.4|) is not well defined, because it may happen that ((U) attains the 
value for a finite value of t and hence that the solution of the Cauchy problem 
(|3.4[) is not a diffeomorphism r : [0, +oo[^ [0, +oo[. 

Conversely, the slow dynamics behave somehow like the component u\ in (|3 . 2(1 . 
namely they are solutions of (|3.ip that satisfy an equation with no singularity in 
it. We refer to Bianchini and Spinolo [4] for a more rigourous definition of fast and 
slow dynamics. 

Before introducing Theorem 13 - 21 we point out that, as a consequence of the 
hypotheses introduced in Section [4l we can perform a local change of coordinates 
such that, in the new coordinates, ((U) is actually one of the components of U, say 
((U) = u\ if U = (mi, . . . UdY- Also, the manifold E — {(ui, ... 0) : u 2 = u 3 = . . .Ud 
is entirely made by equilibria for 

Our result is the following: 

Theorem 3.2. Assume that 

F(6) = o c(o) = o 

and that all the hypotheses introduced in Section [^] are satisfied. Then there exists 
a manifold A4 WS , defined in a small enough neighbourhood of and satisfying the 
following properties. 

(1) The manifold A4 US is locally invariant for (|3.5p and it contains all the 
solutions U(t) decaying exponentially fast to an equilibrium point in 

E = {(ui, 0. . .0) : U2 = it3 = • • -Ud = 0} . 

(2) Any orbit lying on A4 US decomposes as 

U{t) = U slow {t) + U fast {t) + U pert {t). 

Here JJ slow (t) and U* ast (i) are a slow and a fast dynamic respectively, in 
the sense explained before. The perturbation term U pert (t) is due to the 
non linearity and satisfies 

\U pert (t)\ <CC/ /ast (0)C(V(0)) 

for a suitable constant C . 

(3) If U(t) is an orbit lying on M c and Q{U) > at t = 0, then ((U) > 
for every t. Also, the Cauchy problem (|3.4p defines a diffeomorphism 
t : [0, +oo[^ [0, +oo[. 
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Remark 3.3. By direct check one can verify that all the hypotheses introduced 
in Section [4] are verified by (II. 9|) , the equation satisfied by the steady solutions of 
the Navier Stokes written in Eulerian coordinates. 

The application of Theorems 13.11 and 13.21 to (|1.9p is connected to a remark due 
to Frederic Rousset, which is the following. In our hypotheses the Lagrangian and 
the Eulerian coordinates are equivalent. However, steady solutions of the Navier 
Stokes equation written using Lagrangian coordinates are regular since they satisfy 
an ODE with no singularity. As shown by example (|3.2p . in general a solution of a 
singular ODE is not C . If the solutions of (|1.9[) were not C 1 , this would contradict 
the equivalence between Eulerian and Lagrangian coordinates. 

However, Theorems 13.11 and 13.21 tell us that, if we restrict to solutions lying 
on either A4 C or A4 US , then any loss of regularity is ruled out: if v(0) > 0, then 
v(t) > for every t. 

For a different approach to the analysis of the viscous profiles of the Navier 
Stokes equation in Eulerian coordinates see for example Wagner [16] and the ref- 
erences therein. 



4. Hypotheses 

In this section, we introduce the hypotheses we impose on system (|3.5[) . All 
these conditions are satisfied by (|1.9[) , the equation satisfied by the viscous profiles 
of the compressible Navier Stokes equation in one space variable written in Eulerian 
coordinates. 

Concerning the regularity of the maps F and £, we assume that they are C 3 . 
Also, without loss of regularity we can restrict to the case the equilibrium U is 0, 
namely 

F(0) = C(0) = 
We also assume the following conditions: 

(1) The gradient VC(0) ^ 0. 
Let S be the singular set 

S :={U: C(U)=0}. 

Thanks to the implicit function theorem, Hypothesis (1) ensures that in a small 
enough neighbourhood of the set S is actually a manifold of dimension d — 1, 
where d is the dimension of U. 

(2) Let A4° be any center manifold for 

(4.1) % =F W 

around the equilibrium point 0. If \U\ is sufficiently small and U belongs 
to the intersection A4 C D S , then U is an equilibrium for (|4.ip . namely 
F(U) = . 

The reason why we introduce Hypothesis (2) is the following. Consider the linear 
system 

(dui/dr = 112/E 
du 2 /dr = -ui/e 
de/dT = 0. 
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The first component of the solution is 

m(t) = Acos(t/e) + Bsm(t/e), 

where A and B are real parameters. Letting e — > + , we get that in general if t ^ 
there exists no pointwise limit of u\. Note that (|4.2p does not satisfy Hypothesis (2). 
Indeed, the center space is the whole R 3 . However, 

F(U) = 

is not identically zero when e = 0. 

(3) There exists a manifold of equilibria M eq for (14. ip which contains and 
which is transversal to S. 

Let n eq be the dimension of M eq . We recall that the manifolds S and A4 eq are 
transversal if the intersection S PI M eq is a manifold with dimension n eq — 1 (as 
pointed out before, the dimension of S is d — 1). 

(4) For every U e S, 

V<(E/) • F(LQ = 0. 

Hypothesis (4) is necessary if we want to rule out losses of regularity like the one in 
(|2.4[) . Indeed, example (|2.4[) satisfies all the hypotheses stated in this section, but 
Hypothesis (4). 

Thanks to Hypothesis (4) and to the regularity of the functions ( and F, the 
function 

can be extended and defined by continuity on the surface S. 

(5) Let U G S be an equilibrium for (|4.1|) . namely ((U) — and F(U) = 0. 
Then 

G(U) = 0. 
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